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We study theoretically the coexistence of the spin-density-wave (SDW) and superconductivity in electron- 
doped iron-pnictide superconductors based on the two orbital model and Bogoliubov-de Gennes equations. The 
phase diagram is mapped out and the evolution of the Fermi surface as the doping varies is presented. The 
local density of states has also been calculated from low to high doping. We show that the strength of the 
superconducting coherence peak at the positive energy gets enhanced and the one at the negative energy is 
suppressed by the SDW order in the underdoped region. Several features of our results are in good agreement 
with the experiments. 
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The discovery of the iron-based superconducting (SC) ma- 
terials [1] has attracted much attention. These materials share 
some common features of the high-T f cuprates [2], with high 
transition temperature, layered system and similar phase di- 
agram. The parent compounds are with spin-density-wave 
(SDW) or antiferromagnetic (AF) order, and they are not Mott 
insulators but poor 'metals'. By doping either holes or elec- 
trons into the parent compound, the SDW/AF order becomes 
weakened, and superconductivity shows up. Theoretically the 
Fermi surface nesting is proposed to be responsible for the 
presence of the SDW order (3J. Moreover, the mechanism 
for the superconductivity in this system is still unclear while 
many suggested that the correlation effect is still essential and 
the spin fluctuation should be responsible for the supercon- 
ductivity yd. 

A quite intriguing question is how the magnetic (SDW/AF) 
order competes with the superconductivity and whether the 
magnetic order could coexist with the SC order at low dop- 
ing densities. This is still an open issue and the result is 
under debate or depends on the materials. It was reported 
that the magnetic phase and SC phase are totally separated 
without the coexisting region in the fluorine doped material 
CeFeAsOi- A F A according to the inelastic neutron scattering 
(INS) experiments fl5[] . However, it was also reported that the 
magnetic order and the SC order coexist at low doping den- 
sities in the Co-doped material BaFe2- A Co A As2 |6-8]. The 
coexistence of the two orders is also reported in some other 
materials, such as Bai_ y K x Fe2As2 (Hill], Sri_ A K x Fe2As2 ifiotl . 
and SmFeAsOi-jF* 111 111 . Theoretically, it was proposed that 
the coexistence of these two orders is possible and can ex- 
plain some experimental observations 1131 . While relatively 
this subject remains less explored on the theoretical front. 

The motivation of the present work is to fill this void and 
examine this issue for electron-doped samples. We adopt a 
two-orbital model by taking into account two Fe ions per 
unit cell, as proposed in Ref. lfl3tl by one of the present au- 
thors. We map out the phase diagram and show that the 
doping evolution of the Fermi surface based on this model 
is consistent with the angle resolved photoemission spec- 
troscopy (ARPES) experiments on the electron-doped sam- 



ples BaFe2-.1C0jA.s2 iffl [l5ll . We calculate the spatially dis- 
tributed magnetic order and SC order self-consistently based 
on the Bogoliubov-de Gennes (BdG) equations. The SC or- 
der here is chosen to have s ± -wave symmetry which is sup- 
ported by some experimental observations 111611 and theoret- 
ical studies 01, [nil . The obtained SDW phase in which the 
electron spin is ferromagnetic along >'-direction and antifer- 
romagnetic along x-direction is in agreement with the exper- 
iments lfl8l [l9tl . The magnitude of magnetic order decreases 
as the doping increases. The superconductivity occurs as the 
doping 8 > 0.01. The magnetic order coexists with the SC 
order at low doping (0.01 < 5 < 0.12). The local density of 
states (LDOS) is calculated and the signatures of the coexis- 
tence of the SDW and SC orders are also identified. 

We start with an effective model with the hopping elements, 
pairing terms, and on-site interactions, expressed by, 

H = H t +H A + H int . (1) 

The first term is the hopping term, expressed by, 

flf = - 2 (ty> c 5«r C /K r + ~ f ° Yj C l*r c H*r' (2) 

if* jvo~ ifi(r 

where i, j are the site indices and fi, v - 1 , 2 are the orbital 
indices, to is the chemical potential, 
is the pairing term, 

#a= £(Ai^c^c^+/*. C .). (3) 

iftjva- 

Hi nt is the on-site interaction term. Following Ref ll20ll . we 
here include the Coulombic interaction and Hund coupling 
Jh- At the mean-field level, the interaction Hamiltonian can 
be written as ll2lll . 

H inl = U 2^ (n^d-)n itlC7 + U' ^ (n^s-^W 

tViT ■ (4) 
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where « J/JIT are the density operators at the site i and orbital fi. 
U' is taken to be U - 2J H udii . 

Then the Hamiltonian can be diagonalized by solving the 
BdG equations self-consistently, 
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where the Hamiltonian Hi„j V( r is expressed by. 
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The SC order parameter and the local electron density (n^) 
satisfy the following self-consistent conditions, 
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FIG. 1 : (Color online) (a) The intensity plot of the magnetic order at 
zero doping and zero temperature, (b) The magnitude of the SC order 
parameter A and magnetic order m as a function of the temperature 
with the doping density 6 = 0.09. (c) The magnitude of the SC order 
parameter A and magnetic order m as a function of the doping at zero 
temperature, (d) The calculated phase diagram. 



Here Vi^p is the pairing strength and fix) is the Fermi distri- 
bution function. 

The LDOS is expressed by, 

Pi(co) = ^[K^l 2 ^,, " + Ka-MEn + (9) 

where the delta function Six) is taken as F/t:(x 2 + T 2 ), with the 
quasiparticle damping T = 0.004 (the results are qualitatively 
the same for different T). The supercell technique [22] is used 
to calculate the LDOS. 

We use the hopping constant suggested by Ref. lfl3ll . 
namely, 



= h 
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The pairing symmetry is determined by the pairing poten- 
tial Vinjy. We have carried out extensive calculations to search 
for favorable pairing symmetry based on the present band 
structure. Especially, the next-nearest-neighbor intra-orbital 
potential will produce s A 2 v 2-wave (s ± -wave) order parameter, 
i.e., A oc cosfc^cos £ v /(cos£ A +cosfc v ) in the extended/reduced 
Brillouin zone (the Brillouin zone reduces to half of the ex- 
tended one considering two Fe ions per unit cell). Here the 
pairing symmetry we obtained is independent on the initial 
input values and the results is consistent with previous cal- 
culation based on a different two-orbital model ll2lll . In the 
following calculation we will focus on this pairing symmetry. 

Throughout the work, we use the hopping constant ti-4 = 
1,0.4,-2,0.04. fo is determined by the electron filling per 
site n (n = 2 + 5). The on-site Coulombic interaction and 



Hund coupling U and Jh are taken as U = 3.4 and Jh = 
1.3, respectively, consistent with the recent estimation l23ll . 
The pairing potential is chosen as V = 1.2. The numerical 
calculation is performed on 20 x 20 lattice with the periodic 
boundary conditions. An 80x80 supercell is taken to calculate 
the LDOS. 

We plot the spatial distribution of the magnetic order [m, = 
5 ZZfi( n iMl ~ n v4)] at me zero temperature and zero doping 
in Fig. 1(a). As seen, the magnetic spin order is antiferro- 
magnetic along the x direction and ferromagnetic along the 
y direction, corresponding to the (n, 0)l(n,7i) SDW in the ex- 
tended/reduced Brillouin zone. This result is consistent with 
the INS experiments lfl8l [l9ll and previous theoretical calcu- 
lation based on different band structure [21]. There exists an- 
other degenerate SDW state with the spin order is antiferro- 
magnetic along y direction and ferromagnetic along x direc- 
tion. The four-fold symmetry breaking and the presence of the 
SDW order are due to the Fermi surface feature nesting at low 
doping which will be discussed below. The magnetic order 
decreases as the temperature increases or doping increases. 
The amplitudes of the SC order and magnetic order as func- 
tions of the temperature and doping are shown in Figs. 1(b) 
and 1(c), respectively. At the fixed doping, both the magnetic 
order and SC order decrease as the temperature increases and 
two transition temperatures are revealed, as seen in Fig. 1(b). 
The superconductivity occurs at the doping 0.01 and vanishes 
at the doping 0.26, as seen in Fig. 1(c). The magnetic order 
is maximum at zero doping and decreases monotonically as 
the doping increases. The calculated phase diagram is plotted 
in Fig. 1(d). As seen, the magnetic order and SC order co- 
exist in the underdoped region. The magnetic order decreases 
abruptly and a quantum critical point 5 — 0. 12 is revealed. The 
superconductivity appears as the magnetic order is suppressed 
and the SC transition temperature T c reaches the maximum 
as the magnetic order disappears. Our results are reasonably 
consistent with the experiments on the BaFe2- j: Co A: As2 fa-Bal. 
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FIG. 2: (Color online) The doping evolution of the Fermi surface at 
the temperature T = 0.1 (i.e., the superconductivity disappears). The 
(red) bold solid curves and (black) dotted lines in panels (a-c) indi- 
cate the ungapped Fermi surfaces and the magnetic Brillouin zone 
boundary in the magnetic phase, respectively. The (blue) thin solid 
curves in panels (a)-(c) show the normal states Fermi surface by set- 
ting the magnetic order m = 0. 

The real-space mean-field Hamiltonian can be transformed 
to the momentum space because the SC order is uniform and 
the magnetic order has the period 2x1. The Fermi surface, 
defined by the zero energy contours of the quasiparticles, can 
be obtained through the momentum space Hamiltonian. The 
evolution of the normal state Fermi surface with increasing 
doping densities is shown in Fig. 2. 

At zero doping [Fig. 2(a)], the Fermi surface in the normal 
state [(blue) thin lines] contains two hole pockets around point 
r = (0, 0) (labeled as a and /3 bands) and two electron pock- 
ets around point M - (n,7i) (labeled as 5 and y bands). The 
existence of the SDW here represents the pairing of the elec- 
tron and hole. Thus in the SDW phase, the momentum of 
the SDW order should be consistent with the Fermi surface 
nesting momentum connecting the hole pockets and electron 
pockets. As seen in Fig. 2(a), the prime nesting vector (de- 
noted by the dashed arrows) in the reduced Brilliouin zone 
is (n, n). As the doping increases, the hole pockets shrink. 
Then the nesting vector should deviate from (n, n). But with 
the presence of the underline lattice, the nesting vector is still 
pinned at (n, n) and this is what has been observed by experi- 
ments [7], thus the magnetic order will decrease as the doping 
increases. When the doping increases to 5 ~ 0.1 (at T = 0.1), 
there is no nesting vector that can be sustained by the lattice 
and the SDW order disappears. The magnetic Fermi surfaces 
[red (bold) solid curves] as shown in Figs.2(a-c) are repre- 
sented by ungapped pockets along the T — M line or k x = k y 
direction. Parts of the original Fermi surface [(blue) thin solid 
curves] are gapped by the SDW order. The small ungapped 
Fermi surface pockets along the T — M line is consistent with 
experiments on the parent compound BaAs2Fe2 B2411 . The ex- 
istence of ungapped Fermi surface pockets at zero doping [see 
Fig. 2(a)] indicates that the parent compound is in fact not an 
'insulator' a but a poor 'metal'. Here the four- fold symme- 
try is broken in the magnetic phase, due to the presence of 
the magnetic order shown in Fig. 1(a). There exists another 




FIG. 3: (Color online) The LDOS spectra with different doping den- 
sities at the zero temperature. The (red) dash-dotted and (blue) dotted 
lines in panels (b) and (c) are non-selfconsistent results with A = 0, 
and m = 0, respectively. 



degenerate solution with the Fermi surface pockets along the 
k x = -k y direction, corresponding to the (0, 7r)-SDW ground 
state in the extended Brillouin zone. As the doping increases, 
the ungapped pockets become larger and larger and eventually 
cover the whole Fermi surface as the SDW order disappears. 

As the doping increases to 6 — 0.16, and 0.24 [Figs. 2(d) 
and 2(e)], the c-band is filled by electrons completely and the 
Fermi surface of this band disappears, thus only three Fermi 
surface pockets are left. This feature is consistent with the 
experiments 11411 . As the doping density increases to about 
0.30 [Fig. 2(f)], the Fermi surface pocket of yS-band will also 
disappear and only two electron-like Fermi surface pockets 
are left, which is also consistent with the experiments 1151 . 
On the other hand, the electron pockets expand a little as the 
doping increases but relatively the electron pockets depend 
weakly on the doping, also in agreement with the ARPES ex- 
periments 11141 Il5ll . 

We also calculate the LDOS spectra according to Eq. (9). 
The spectra at different doping densities are shown in Fig. 3. 
At zero doping [Fig. 3(a)], the density of states shows four 
well-defined coherence peaks due to the SDW order. The 
maximum dip of the spectrum is at the chemical potential or 
zero energy. Note that there are two coherence peaks at nega- 
tive energies and two at positive energies. The splitting of the 
coherence peaks is caused by the inter-orbital coupling t\. As 
the doping increases [Figs. 3(b) and 3(c)], the maximum dip 
of the LDOS due to the (n, n) SDW order is nearly indepen- 
dent of doping and always pinned at the chemical potential 
of zero doping. This is similar to the case of d-density wave 
in cuprate superconductors 12511 . As we define the zero en- 
ergy at the chemical potential, the maximum dip of the SDW 
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spectrum will shift to the left or the negative energy as dop- 
ing increases. While in the SC state, two SC coherence peaks 
show up, and the mid-gap point is always located at the zero 
energy or the chemical potential of the electron-doped system. 

We now discuss the LDOS spectra in the coexisting region. 
In this region, the pure SC state and the pure SDW state have 
higher free energy. These two pure states can be obtained by 
setting the magnetic order m = and the SC order is cal- 
culated self-consistently and setting the SC order A = and 
the magnetic order is obtained self-consistently, respectively. 
The results of the pure SC and SDW states are represented 
by (blue) dotted and (red) dash-dotted lines, respectively in 
Figs. 3(b) and 3(c). In the pure SDW state, four SDW peaks 
are obtained with the maximum dip at the negative energy. In 
the pure SC state, two SC coherence peaks are seen with the 
mid point of the gap at the zero energy. In the coexisting state, 
the structure of the spectra inside the SC gap is practically the 
same to the spectra of the pure SC state. As shown in Figs. 3(b) 
and 3(c), we can see clearly the previous SDW peaks outside 
the SC gap. But the SDW peaks inside the SC gap in the pure 
magnetic state cannot be seen in the coexisting state. The ad- 
ditional peaks outside the SC peaks, positioned by the (red) 
vertical arrows can be regarded as one signature of the coex- 
istence of the magnetic and SC orders in electron-doped ma- 
terials. As doping increases, only negative SDW peaks would 
appear, as seen in Figs. 3(c) and 3(d). It is needed to point out 
that such structures outside the SC gap due to the SDW order 
as shown in Figs.3(a-d) so far have not been clearly identified 
by the experiments. 

Another significant feature caused by the magnetic order 
can be seen from the intensities of the SC coherence peaks 
[see Figs. 3(b) and 3(c)]. In pure SC state, the intensity of the 
SC peak at the negative energy is higher than that at the posi- 
tive energy for all the doping densities we considered. These 



SC peaks are positioned by the parallel (blue) dotted arrows. 
In the coexisting region, as we mentioned above, the SDW 
spectra shift to the negative energy so that the SC peak at the 
negative energy is within the SDW gap and that at the positive 
energy is outside the SDW gap. Thus the intensity of the SC 
peak at the positive energy is enhanced and the one at the neg- 
ative energy gets suppressed by the SDW order. As a result, 
in the coexisting state, the intensity of the SC coherence peak 
[denoted by the (black) solid arrows] at the negative energy is 
lower than that at the positive energy. The asymmetry disap- 
pears near the optimal doping (« 0.105), as seen in Fig. 3(d). 
In the overdoped region [see Figs. 3(e) and 3(f)], the asym- 
metry occurs again but with the intensity of the SC peak at 
negative energy becoming stronger. This feature has recently 
been confirmed by the scanning tunneling microscopy (STM) 
experiments on BaFe2-j;Co v As2 B26I1 . 

In summary, we have examined theoretically the coexis- 
tence of the SDW and SC orders in electron-doped iron- 
pnictide superconductors based on the two orbital model and 
BdG equations. The phase diagram is mapped out and the co- 
existence of the SDW and SC orders occurs at low doping. 
The evolution of the Fermi surface as the doping varies is pre- 
sented and the results agree with several ARPES experiments. 
The LDOS has also been calculated from low to high doping. 
The signatures of the SDW order are identified and is consis- 
tent with the recent STM experiment. It is important to point 
out that the present results are critically dependent on the two 
orbital model suggested in Ref. Ill 311 . 
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